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Grover's algorithm for quantum search can also be applied to classical energy transfer. The 
procedure takes a system in which the total energy is equally distributed among N subsystems 
and transfers most of the it to one marked subsystem. We show that in a relativistic setting the 
efficiency of this procedure can be improved. We will consider the transfer of relativistic kinetic 
energy in a series of elastic collisions. In this case, the number of steps of the energy transfer 
procedure approaches 1 as the initial velocities of the objects become closer to the speed of light. 
This is a consequence of introducing non-linearities in the procedure. However, the maximum 

■ attainable transfer will depend on the particular combination of speed and number of objects. In 
the procedure, we will use N elements, like in the classical case, instead of the log2(N) states of the 

, quantum algorithm. 

O: 

'. INTRODUCTION 

Physics gives new ways to analyse information transmission and computation. The traditional theories of informa- 
Oh| tion and communication implicitly assume an underlying classical world, with no regard to modern physical theories. 

■ However, the recent physical analysis of information has considerably enriched both fields. Taking quantum mechanics 
into account in communications has led to quantum cryptography [l| 0] an d the applications to information processing 
include efficient algorithms for factoring [3( among others [J]. 

q-i But, while quantum information has become an active field of research, there are few results that incorporate the 
l— ~~ '■ principles of the theory of relativity into communications and information processing. There are already promising 
I | studies that include relativistic constraints into new communication protocols a, |6j. In particular, relativistic bit 
commitment protocols offer results with no classical or quantum counterpart 0, 8|]. In spite of these advances, there 
| has been no significant work on the possibility of relativistic computing. 
\D . In this paper, we take a variation of Grover's algorithm for quantum search and translate it into a relativistic 
setting. This variation can be though of as an effective procedure for energy transfer between N classical objects. 
In the classical case, the algorithm needs a number of steps proportional to y/~N to transfer the energy to a marked 
object. We show that special relativity predicts a smaller number of steps, which can be reduced even to a single step 
when the speed of light, c, is approached. 



On 
O 



ON 

o 



- I - 1 

>< 



C3 



GROVER'S ALGORITHM 



The basic form of Grover's algorithm for quantum search gives a procedure to stand out a marked element out of 



N possibilities [9( . The algorithm relies on an oracle that can detect and invert the phase of a particular target state 
\t). The starting point is a uniform superposition \s) of all the possible states. With a sequence of oracle calls and 
diffusion operators that redistribute the probability amplitudes of a state among the others, we can gradually transfer 
the probability amplitude from the non-marked states to the target. After a number of iterations of the order of yN, 
there is a high probability of finding the desired state in a measurement. 

The evolution can be more easily understood if we define a state \t±) that collects all the non-marked states from 
the original superposition. This state is orthogonal to \t) and close to |s) in the original Hilbert space. If we observe 
the evolution of |s) in the two-dimensional space spanned by \t) and \t±), Grover's algorithm reduces to a series of 
rotations that take the state from |s) to \t) [l0(. 

A similar procedure can be defined to transfer energy between classical oscillators [ll[. In the initial state, the 
total energy of the system is equally distributed among N oscillators. During the Grover procedure, the energy of a 
marked oscillator will go from being of the order of ^th of the total energy to being of the order of the total energy. 
Now, the diffusion operator redistributes energy instead of probability amplitude. The crucial difference is that in the 
classical case we can no longer use a quantum superposition. This results in an increase in the space complexity of the 
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FIG. 1: An iteration of the Grover energy transfer. In the initial state, the big ball moves to the right with velocity vi and 
the small one with velocity vi. First, the small ball rebounds from the wall in an operation equivalent to the oracle sign shift 
of quantum search. After that, both balls collide. It is in the collision where the energy is transferred between the balls. This 
step is equivalent to the diffusion operation of Grover search. As a result of the collision, the balls move with velocities v[ and 
v' 2 . This final state is the initial state of next iteration. 

procedure. In the quantum case we could encode N states in log2N two-level systems. Now, we need N oscillators. 
Apart from that, the time complexity of the procedure (the number of steps) is the same in both scenarios. 

Grover energy transfer can also be better studied by defining two different modes. The first one corresponds to the 
marked mode to which we want to transfer the energy of the global system of N oscillators. The second mode is the 
collective mode of all the non-marked oscillators. 

We will discuss a simplified case based on elastic collisions [l2| • The basic model will be the same for both classical 
and relativistic Grover energy transfer. We will consider N bodies of normalized mass 1 which, to simplify the study, 
will be grouped in two balls: a big ball of normalized mass M = N — 1 and a small ball of normalized mass 1. Both 
balls are moving to the right with the same initial speed Vq. Our algorithm will be designed to transfer all the kinetic 
energy of the composite system to the small ball. The transfer will occur in a series of elastic collisions. 

The situation is summed up in Figure [TJ We begin with both balls moving to the right (with positive speeds vq). To 
simplify things, we take the smaller ball to be a little ahead of the big one. If needed, this separation could be part of 
the oracle phase. The small ball will meet a rigid wall that will act as the oracle in Grover's algorithm. Now, instead 
of a 7r phase shift, we have a change of sign in the velocity when the small ball rebounds. After the rebound, the 
balls will collide. Collisions play the role of the diffusion operator. They redistribute the energy into all the modes, 
in this case, the kinetic energy of the small and big balls. This will be the setting of both the classical and relativistic 
Grover energy transfer procedures. The small ball will start with one Nth of the total kinetic energy and, at the end 
of the procedure, its fraction of the total kinetic energy will be close to one. We suppose that during the whole series 
of collisions, the big ball does not meet the wall and still moves to the right. This is not strictly necessary, but makes 
the analysis more direct. 

In the classical case, the evolution of the system after each collision is derived from conservation of energy and 
momentum. This evolution, when written in a matrix form, corresponds exactly to the evolution operator of quantum 
search and requires the same number of steps [lij ]. Notably, this does not depend on the initial velocity. 

Here, we generalize the situation to relativistic speeds. Our analysis is based on special relativity and will pay no 
regard to any possible gravitational effects. All the given parameters are considered in the frame in which the wall is 
stationary. 

We assume a series of elastic collisions. The evolution can be deduced from the conservation of four-momentum. 
We consider a simplified case in which there is only linear movement in the x direction. All the velocities given will 
be taken to be in the x direction and are normalized to the speed of light, using Vi as a shorthand for In this case, 
conservation of four-momentum can be expressed with only two equations (for conservation of relativistic energy and 
^-momentum) . Additionally, as we are dealing with elastic collisions, the rest mass of each ball will be conserved. 
The evolution after a collision is given by 



for the conservation of energy and x-momentum, respectively. In the equations, 7$ represents the Lorentz factor 
, 1 , corresponding to velocity The left hand side of the equations define the four- momentum components 

before the collision and the right hand side the situation after the collision. Variables with subindices 1 refer to 
parameters of the big ball and those with subindices 2 correspond to the small ball. 

The resulting non-linear system of equations can be more easily solved by adding and subtracting both equations 



Mji + 72 = + 7 2 

M71U1 + 7 2 «2 = Mf[v'i + 1W2 



(1) 

(2) 
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to give the equivalent system, 



Mq! + q 2 = Ma[ + a' 2 (3) 
Ml Ml 

-T + -T, (4) 



0t\ 0i2 



where a, = *J ^z^f. ■ Now, we only need to solve the quadratic equation Aa 2 + Ba' 2 + C — 0, with A = + 

B = — [M( ^ + ^ )+2] and C = Mai + «2- From the two possible solutions, we take the one that corresponds to 
the big ball still moving to the right and the small ball rebounding from it. With this value of a' 2 , we can recover 

a 2 1 

a'i = oti + jj{a>2 — Q^)< The final velocities can then be calculated from u,; = "i, x . 

While we could study the energy transfer using velocity, we can understand the evolution better by looking into the 
relativistic kinetic energies of the balls. We have total energies E\ = M71 and E 2 — 72, but, as masses are conserved, 
the only energy that can be transferred is kinetic energy K\ = M71 — M and K 2 = 72 — 1. We are not interested in 
the particular values of these energies, but in their relationship to the total kinetic energy Kt- Kt is conserved all 
throught the energy transfer procedure and can be deduced from the initial condition Kt = K[ m + K 2 m — TV 70 — N. 



Results 



Reduction in the number of steps 

Figure [2] shows the number of steps needed for the maximum energy transfer. As in Grover's algorithm, if we do 
not stop the procedure, the kinetic energy of the system will alternatively go from the big to the small ball. We have 
only considered the first maximum of kinetic energy of the small ball (the fastest energy transfer). In the initial state, 
jp- is exactly -j^. When we stop the procedure, K 2 is of the order of Kt- The results are given for initial velocities 
vq = 0.001,0.01,0.05,0.1,0.3, and 0.8. The number of steps is plotted along the asymptotic limit of steps for Grover 
search, f 0. 

The graph shows some representative cases of relativistic Grover energy transfer. The most significant conclusion 
is that, for this relativistic setting, the procedure is more efficient than in the classical case. The reduction in the 
number of steps is greater for higher velocities. As expected, for small velocities, the deviation from the classical case 
is small. However, as we come closer to the speed of light, the number of steps is dramatically reduced. For a high 
enough initial velocity vq, the maximum transfer can occur even in a single collision. 

We can estimate the initial velocity needed for this single step transfer from four-momentum conservation. In our 
only collision, equations (Q]) and © become 

N l0 = M + 7 2 (5) 
N~f Q v -2j v = j' 2 v 2 . (6) 

The left hand side shows initial energy and momentum and the right hand side corresponds to a perfect transfer in 
which the big ball stops (v[ = 0) and all the kinetic energy goes to the small ball. If we take vq = 1 — x and v 2 = 1 — y, 
for x,y « 1 and assume N >> 1, we can subtract the equations and find a single step velocity Vq s — 1 — ™. 
Numerical simulations show maximum energy transfer in a single collision for a given TV and values of vq around Vq s . 

Even for very small initial velocities, there will be a value of N for which relativistic effects will become relevant. 
This can be best noticed looking at the graph of Vq = 0.001 in Fig. [21 For moderate values of M the deviation from 
classical is negligible. However, for high values of M the curve starts to separate from the classical prediction. 

We can explain the separation looking into the Taylor expansion of K\ for small initial velocities, K\ = \Mvq + 
|Mwq + 0(vq). In most problems with a small vq, we could neglect all the terms but the classical kinetic energy. 
However, for high values of M, the correction term can be comparable to the other relevant energy of the problem, 
K 2 . We have taken K\ for the initial state, when v\ = vq, as a limit. During Grover energy transfer, v\ becomes 
smaller at every step as the big mass losses energy. The initial state gives us the maximum deviation of K\ from the 
classical approximation. 

For a given problem (a particular value of N), we can define a breakpoint initial velocity Vq at the value at which 
the first correction term of K\ is equal to the classical approximation of K 2l ^Mv^ b — \v^ b . This happens when 
Vq = J^ M . Velocities much below Vq will follow the classical Grover energy transfer, while velocities well above «o6 will 

show a clear deviation. Similarly, for a given vq, we can define a breakpoint M b = This estimation corresponds 
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FIG. 2: Number of iterations before the maximum transfer of energy. As the initial speed increases, the number of collisions 
needed to transfer most of the kinetic energy to the smaller ball is reduced. 

well to the behaviour shown in Fig. [2] and is best noticed for Vq = 0.001 (with iV& » 1.3 • 10 6 ) and vq = 0.01 (with 
N b w 1.3- 10 4 ). 



Transferred fraction of energy 

We can also look at the fraction of the total energy that is transferred at the end of the procedure. Figure [3] shows 
that, while in the final step the energy of the small ball is still in the order of the total energy, there can be a significant 
variation for different initial conditions. 

This behaviour is not exclusive of the relativistic case. In the quantum and classical Grover algorithms there is also 
a point at which an additional collision cannot increase the energy of the marked element any more. In the Grover 
procedure, energy is transferred in discrete steps. The final energy of the small ball comes from adding these steps 
and cannot exceed the total energy. Steps with an energy transfer above the difference Kt — K2 are only allowed if 
they take energy from the small to the big ball. 

This can be compared to the end of a Game of the Goose game. If the die roll exceeds the exact value needed to 
reach home, the piece rebounds and moves back the excess number of positions. 

In a Grover procedure, a total or almost total transfer can be achieved if the discrete transfer steps carry a low 
energy and we have a greater precision. When each step transfers a higher energy /probability, as it happens in both 
classical and quantum Grover algorithms for small values of N, the initial value of N will determine the maximum 
final amplitued. The difference between the maximum transfer values vanishes as N grows and the steps are more 
precise. For the same reason, the varying maximum amplitude effect becomes more marked at relativistic speeds. 
When the initial velocity is high, the energy transferred at each step is greater and we have a rougher approximation 
to an exact transfer of all of the kinetic energy. 

This means that choosing an appropriate pair of values vq and N can improve the efficiency. There will be some 
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FIG. 3: Fraction of the total kinetic energy Kt that the small ball has at the end of the Grover energy transfer. 



TABLE I: Comparison of the efficiency of equivalent quantum, classical and relativistic Grover procedures. 





Number of elements 


Number of steps 


Quantum 


log 2 (N) 


o(Vn) 


Classical 


N 


o (Vn) 


Relativistic 


N 


l(v -» 1) 



initial values for which we can obtain a perfect fit to the maximum energy (an exact die roll). 

ANALYSIS AND FUTURE WORK 

We have presented results showing that the effective procedure for Grover energy transfer can be more efficient 
when taken from a classical to a relativistic setting. Table U gives a comparison between the time and space resources 
required by the different Grover procedures in the quantum, classical and relativistic cases. The space resources 
(number of elements) involved in the classical and relativistic cases are both of the order of N. In the quantum 
Grover procedure, where we transfer probability amplitude, quantum superposition allows for a greater compactness 
as log 2 N two-level states are enough to encode N values. 

The time complexity (number of steps) grows as y/W in both the quantum and the classical procedures. This bound 
has been shown to be optimal However, the nonlinearities that appear in the relativistic version of the procedure 
alter the behaviour and allow a faster transfer. For velocities close to the speed of light the transfer can even be done 
in a single step. 

These results show that there are physical algorithmic procedures that are more efficient in situations where rela- 
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tivistic effects apply. This suggests a new field of Relativistic Computation, parallel to Quantum Computation. In 
this particular case, direct application is doubtful, although the Grover energy transfer procedure could be of interest 
in particle colliders, as long as elastic collisions could be achieved. 

Nevertheless, the procedure could be implemented with other nonlinear analogues. There have already benn analog 



simulators of interesting relativistic effects like Hawking radiation in optical fibre systems LjJ . Similarly, relativistic 
Grover energy transfer could be simulated with optical nonlinear systems inside optical fibre or microring resonators. 
This also opens the door to a new line of research on nonlinear Grover procedures in different systems like, for instance, 
a series of electronical oscillators with nonlinear coupling. 

There are also many possible extensions to the presented Grover transfer procedure. The most obvious addition 
is checking the procedure's efficiency when we have more than one marked object. In a quantum Grover algorithm 
with k targets, the number of steps can be divided by a factor of y/k [l5j |. We have already performed preliminary 
simulations for more than one marked object and, while there is an improvement for high velocities when compared 
to the classical transfer, the advantage reduces as the number of marked objects increases. A detailed analysis of this 
data will be presented elsewhere. 

Additionally, it would be desirable to have at least a rough approximation to the number of steps the procedure 
requires for each initial velocity. One of the greatest challenges in dealing with relativistic Grover energy transfer is 
finding intuitive explanations for the results. The nonlinearity of the equations obscures an easy interpretation. Our 
work in progress includes a search for closed expressions for the assymptotic limit of steps and combinations of vq and 
N that optimise the transfer. Curiously, the classical analogue of Grover algorithm gives a simple physical argument 
of why the optimal number of steps should be of the order of \/N [lfj . 

Relatedly, it is still to be proven if the relativistic procedure is also optimal. The fact that we recover the \ZjV bound 
the at the classical limit (vq — > 0) suggests that it is, but there could be a different behaviour for higher velocities. 

Finally, in a more general programme for Relativistic Information, it appears the exciting possibility of adapting, 
or even creating, new different algorithms considering both special and general relativity. 
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